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Abstract 

We study the leptogenesis scenario in models with multi-Higgs doublets. It is pointed out 
that the washing-out process through the effective dimension five interactions, which has 
not been taken into account seriously in the conventional scenario, can be effective, and 
the resultant baryon asymmetry can be exponentially suppressed. This fact implies new 
possible scenario where the observed baryon asymmetry is the remnant of the washed 
out lepton asymmetry which was originally much larger than the one in the conventional 
scenario. Our new scenario is applicable to some neutrino mass matrix models which 
predict too large CP-violating parameter and makes them viable through the washing- 
out process. 
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The origin of the observed baryon asymmetry, the ratio of the baryon density to the 
entropy density 

Fb = ^1 = ( 0. 6 -1)x10- 10 , 

s 

is one of the major problems in cosmology. Various scenarios of baryogenesis have been 
discussed and leptogenesis fU [H is one of the attractive scenarios. Providing U(1)b-l 
breaking in the original theory, lepton number is generated through the CP-violating 
out-of-equilibrium decay of the right-handed Majorana neutrinos. A part of this lepton 
number is converted into the baryon number via the (B+L) violating sphaleron process 
being in thermal equilibrium Q. As a result, the baryon asymmetry in the universe is 
generated. 

The leptogenesis is also an interesting scenario in the following point of view. Amount 
of the generated lepton asymmetry (baryon asymmetry) is related to the neutrino Dirac 
mass matrix, and thus related to neutrino oscillation data || through the see-saw mech- 
anism |J which naturally explains the smallness of the neutrino masses. For detailed 
discussion, some concrete models are necessary. There have been many works on the 
models of neutrino mass matrix, where the relationship between the neutrino oscillation 
data and the observed baryon asymmetry are discussed in detail. 

While the conventional leptogenesis scenario has been discussed for one Higgs doublet 
model, in this letter, we study the leptogenesis in the models with (up-type) multi-Higgs 
doublets. In multi-Higgs models, one-to-one correspondence between the Yukawa cou- 
plings and the neutrino Dirac mass matrix is no longer valid. This fact leads to new 
possible scenario of the leptogenesis, where the washing-out process through the dimen- 
sion five interactions plays the crucial role for the resultant baryon asymmetry. Since the 
effect of this washing-out process has not been considered seriously, the scenario is worth 
investigating. 

Let us first briefly review the conventional leptogenesis scenario. In the following, 
our discussion is always based on the effective Lagrangian at energies lower than the 
right-handed neutrino masses such that 

£n = -//, ; 77~o.\; ; - -Y^NfMiNi + h.c. , 

i 

where i,j = 1, 2, 3 denote the generation indeces, h is the Yukawa coupling, l L and are 
the lepton and the Higgs doublets, respectively, and Mj is the lepton-number-violating 
mass term of the right-handed neutrino iVj (we are working on the basis of the right- 
handed neutrino mass eigenstates). For simplicity, we assume the hierarchy among the 
right-handed neutrino masses, Mi <C M2 <C M3, in the following. 

The lepton asymmetry in the universe is generated by CP-violating out-of-equilibrium 
decay of the heavy neutrinos, N — ► li<p and N — > l^tp. The leading contribution is given 
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by the interference between the tree level and one-loop level decay amplitudes, and the 
CP-violating parameter is found to be J7j[] 



\— -£> \(tfh%] {f(Mf/M 2 1 )+2g(Mf/M!)} 



Here f{x) and g(x) correspond to the vertex and the wave function corrections, 

'1 + x s 



f(x) = yfx 

g(x) i 



1 - (1 +x)\n 

X 



X 



2(1 -x) 



(2) 



respectively, and both are reduced to ~ — 2\/x ^ ^" have assumed that the 
lightest Ni decay dominantly contributes to the resultant lepton asymmetry. In fact, 
this is confirmed by numerical analysis in the case of hierarchical right-handed neutrino 
]. Using the above e, the generated Yb is described as 



masses 



Y, 



B 



~ —d 

9* 



(3) 



where g* ~ 100 is the effective degrees of freedom in the universe at T ~ Mi, and d < 1 
is so-called the dilution factor. This factor parameterizes how the naively expected value 
Yb ~ e/ g* is reduced by washing-out processes. 

We can classify the washing-out processes into two cases with and without the external 
leg of the heavy right-handed neutrinos, respectively. The former includes the inverse- 
decay process and the lepton-number-violating scatterings mediated by the Higgs boson 
such as N + Il ^qji + qL, where qi and qn are quark doublet and singlet, respectively. 
The latter case is the one induced by the effective dimension five interaction, 



h 2 



(4) 



after integrating out the heavy right-handed neutrinos. Here we have described the in- 
teraction symbolically omitting the generation indeces. This term is nothing but the one 
providing the see-saw mechanism 0. The importance of this interaction was discussed 
in |lUfl , where the interaction was shown to be necessary to avoid the false generation of 
the lepton asymmetry in thermal equilibrium. While numerical calculations |J are 
necessary in order to evaluate the dilution factor precisely, Y B ~ e/g* roughly gives a 
correct answer, and the washing-out process is mostly not so effective. Note that this is 
the consequence from the current neutrino oscillation data as explained in the following. 

1 Through out this letter, our notations are all based on Ref. 0. 
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The condition for a washing-out process to be effective is roughly given by 

T(T ~ Mi) > H(T ~ Mi) - , (5) 

where F denotes the decay width or the thermal-averaged cross section times number 
density for the lepton-number-violating scatterings, H is the Hubble parameter, and 
Mp ~ 10 19 GeV is the Planck mass. This condition leads to the lower bound on the 
light Majorana neutrino mass eigenvalues.^ For simplicity, let us assume (h)h)uv 2 / Mj_ = 
{m} D mp,)n/ Mi ~ m ui through the see-saw mechanism, where v = 174 GeV is the vacuum 
expectation value (VEV) of the Higgs doublet, and m ui is the mass eigenvalue of the light 
Majorana neutrino in the z-th generation. 

Considering that the (inverse) decay width is given by 

r D ~ c (h ] h) ll M l ~ c m vl Ml/v 2 , (6) 

where c ~ 10~ 2 is the phase space factor, we can find m v \ > O(10 _3 eV). Although 
estimation for the Higgs mediated scatterings is much more complicated, we obtain the 
similar result. To evaluate the condition for the washing-out process due to the dimension 
five interaction needs knowledge of the flavor structure in the neutrino sector. For sim- 
plicity, assume that the dominant contributions are coming from each single iVj exchange 
for i fixed. In this case, we find the cross section 

^cE(W~fiE^A (7) 

i i 

and the thermal-averaged cross section times number density is given by T(T ~ Mi) ~ 
aMf. This leads to the bound such as > 0(eV 2 ) ( 1Ql " A ^ eV ). Note that this re- 

sult is incompatible with the neutrino oscillation data []5] in the case of hierarchical light 
neutrino masses, since Y,i m li ~ ^ m e ~ 10~ 3 eV 2 , where Am^ is the neutrino oscillation 
parameter relevant for the atmospheric neutrino anomaly. For the degenerate case it may 
be allowed, but may have a conflict with the constraint J2i m vi < 1.8 eV due to observa- 
tions of the large scale structures in the present universe ||11|| . The possibility Mi ^> 10 10 
GeV would be unlikely in the view point of the re-heating temperature of the inflationary 
universe. The condition for the inverse-decay can be consistent with the neutrino oscil- 
lation data, and can be effective. In this case, there is the useful approximation formula 
fT2|l : d ~ l/{k (logfc) ' 6 ) for k > 1, where k ~ T p>/ H(M\). Considering the relation 
d ~ 1/k oc l/m v \, the neutrino oscillation data and the cosmological constraint, we can 
expect d > 10~ 3 . As a result, both of the washing-out processes are not so important for 
the resultant baryon asymmetry in the conventional leptogenesis scenario. 

Now let us study how the above conclusion is changed in models with the multi-Higgs 
doublets. For simplicity, we consider the model with two up-type Higgs doublets, iff 

2 See Eqs. ( p4| ) and (|2^) for more precise evaluation. 



3 



and H% . This simple model includes all the essential points of our discussion. The Dirac 
neutrino mass matrix is written by the sum of two Higgs doublets Hf, 

m D = h^H?) + h 2 (H%) = (h x cos/? + h 2 sin/?) v , (8) 

where v and (3 are defined as v 2 = (H±) 2 + (H%) 2 and tan/3 = [H^ )/ {Hf) , respectively. 
In the following analysis, we take v — 174 GeV as in the conventional scenario, for 
definiteness. It is straightforward to extend all the formulas necessary in the leptogenesis 
scenario to the ones in the multi-Higgs case. For example, the CP-violating parameter 
Eq. ([!]) is modified to be the sum of 

e a = - ri- - EE Im Mh^MKyj+^hiK^Mhh)^ ¥f (9) 



11 b j 



Mi 



for Mi <C Mj, where the first (second) term corresponds to the vertex (wave function) 
correction. The crucial modification in the multi-Higgs models is that the direct rela- 
tion between the Yukawa couplings and the Dirac mass matrix to be compared with 
the experimental data is no longer valid, because of the new parameter (3 introduced. 
Even if the mass matrices, h\{H^) and /^(.ffg), are fixed by some predictive model 
and/or experiments, there is still freedom to change the Yukawa couplings according to (3. 
Since general discussion about the multi-Higgs case is very complicated, we assume that 
hi(Hi) <C ^2(^2) an d hi -C ti2 to make the essential point of our discussion clear. These 
assumption reduces our formulas to the one in the conventional scenario with identifica- 
tion h = ti2- However, note that the quantity defined as TTTd = hv is not the Dirac mass 
to be compared with the experiments, and no longer suffer from the neutrino oscillation 
data. To avoid confusion, we use the notation, {h)h)av 2 / Mi = (m^ D m D )ii/ Mj ~ m U i, in 
the following. Note again that m Vi is not the physical mass eigenvalue. 

In the conventional scenario, the direct correspondence between the Yukawa coupling 
and the Dirac mass matrix allows us to describe the CP-violating parameter as 

E Im [(^^)y]^- (10) 

i 1 



167r(mV?2 D )ii 

by using the Dirac mass matrix. Suppose that is fixed by some mass matrix models 
and/or experiments, we can evaluate e. The resultant e and Y# ~ e/g* may reveal a large 
discrepancy against the observed baryon asymmetry. However, note that, in our case, the 
above formula is modified as e -> e/sin 2 /? through nip — > = m^/sin/?. Therefore, 
in the case with too small e, we can easily reproduce the observed baryon asymmetry by 
adjusting small sin/3 appropriately. This is the straightforward consequence according 
to the introduction of new free parameter (3. On the other hand, if e is too large, there 
seems to be no hope to obtain the baryon asymmetry consistent with the observation. 
However, we will show that, in this case, the washing-out processes can play an important 
role so as to reduce the resultant Y B to the observed values. Since we no longer suffer 



4 



from the neutrino oscillation data, the Yukawa coupling h = h 2 can be taken to be large 
enough for the washing-out process to become very effective. This is a new scenario in 
the leptogenesis which has neither been considered seriously nor be able to be realized in 
the conventional scenario because of the constraints from the neutrino oscillation data. 

Now let us discuss the washing-out processes in detail by integrating the Boltzmann 
equations out. In the conventional scenario, numerical analysis have been done in detail 
||, and all the quantities we need in this analysis are collected there. Thus, we do not 
discuss about derivations of formulas etc. in this letter. All our formulas are based on 
Ref. ||. We are interested in the parameter region where the washing-out processes are 
effective. Although the conventional analysis has been performed also on this parameter 
region, it has not been seriously taken into account, because of the incompatibility with 
the neutrino oscillation data as discussed above. 

The Boltzmann equations are written by |§ 



dY Nl 
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(Y Nl 


- 1 
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sH(M x ) 


l yeq 
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2 Y f rt ,„ (7JV + iNt) 
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where z = M±/T, and s is the entropy density. Here 7d is given by 

eg K X (z) 



D 



with the usual decay width of N\ in the rest system, 



D 



m vl Mj_ 

,2 



(11) 
(12) 

(13) 
(14) 



87T 87T f z 

the modified Bessel functions K iy and the number density = Mf /(n 2 z)K2(z). Here 
7jv (jm) corresponds to the lepton-number-violating scattering process l(f> <-» l(f> (11 <-> 
<fxp) through the effective dimension five interaction of Eq. (H). We have omitted the 
Higgs mediated processes for simplicity. These processes have an effect similar to that of 
the inverse-decay, since both are classified in the same washing-out process with the N\ 
external leg and have amplitudes proportional to (Jvh)n. For two body scattering, 7 is 
given by 



T 



dsa(s)yfsKi 




(15) 



647T 4 J(m a +m b )2 

with the reduced cross section a(s). Assuming, again, that the dominant contribution is 
due to the single iVj exchange for fixed i, the sum of the reduced cross section for the 
scattering l(p <-> l<p and 11 (jxf) is given by 



a N+Nt {s) = cr N (s) + a Nt {s) 



m 



v3 



2tt 



(16) 
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Tam-w* = In + 1m ~ — " 3 4 1 z 6 . (17) 



Here, we have further assumed the hierarchy rn v \^ <C m u3 for simplicity. This also means 
that the washing-out process mediated by N 3 is dominant. Substituting the reduced cross 
section into Eq. (0), we obtain 

4-7T V 4 

This result is a good approximation for the region Mz < T < Mi <C M3 which we are 
interested in, where Mz is the Z-boson mass. 

Now we are ready to numerically analyze the Boltzmann equations. Given e, m u i and 
m U 3, the Boltzmann equations can be easily integrated out. For very small rn v \ and m v3 , 
we can find Yb-l ~ e l 9*- Our interest is on the parameter region where the washing-out 
processes play the important roles. Now we are free from the neutrino oscillation data, 
and can take large input values for m v \ and m v3 . Let us consider two spacial cases. One 
is that the washing-out process through the inverse-decay is dominant. This case has 
been taken into account in the previous works with the range of the resultant dilution 
factor d > 1CT 3 . The other is that the resultant Yb-l is dramatically reduced by the 
washing-out through the dimension five interactions. For each cases, Yb-l is depicted in 
Fig. 1 and Fig. 2, respectively, as a function of Log lo (z).0 Here, we have taken the input 
parameter e = 10~ 3 as an example. The resultant dilution factor d is independent of this 
input in this numerical analysis. We can see that there is the dramatic reduction of the 
resultant baryon number in each graph. 

There are useful approximation formulas for two cases discussed above. The Boltz- 



mann equation of Eq. (|12|) can be rewritten as 

dY B - L 



dz 

where 



+ P(z)Y B - L = Q(z) , (11 



P{Z) = sH(M )Y eq {h D + 2 ^ N+m ) = Pd ^ + Pn ^ ( 19 ) 



sH(M 1JU 
with 



- s£$fa* K M' (20) 

3ml 3 Mf , 

Pn{z) ~ 4tt3#(M 1 K* ' (21) 



and 



Yi 



Q ^Z^mrA^r 1 ) lD ' (22) 

3 The region of the Yukawa couplings in the first case is limited. For large Yukawa couplings, the 
washing-out process through the dimension five interactions usually becomes important, and dominates 
(see Eqs.(E7p and (Eq), and compare them). 
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Here we have used sYf q = -\M 3 z 3 . The solution of the above Boltzmann equation with 
the initial condition Yb-l(O) = is given by 

Y B -l{z)= l Z dxQ{x)e-^ p{y)dy (23) 
jo 

with the solution of Eq. (fill) . We can find that Yb-l(oo) ~ tj 1 g* if P(z) <C 1. 

Noting that the out-of-equilibrium decay begins roughly at T ~ M\, the condition 
P(l) > 1 gives a reasonable condition for the washing-out processes to be effective. For 
the above two cases, we can obtain the conditions on the light neutrino mass eigenvalues 
such that 

Pd(1) > 1 -> m rf >4x l(T 3 eV , (24) 



10 10 GeV\ 



1/2 



Pn(1) > 1 -> m, 3 > 2 eV ^ — j . (25) 

These results are in good agreement with the numerical calculations as can be seen in 
Fig 1 and 2. 

The washing-out processes play the important roles in the case P(l) 1. Using the 
method of the steepest decent in Eq. (p3|) , we can find 



xr / \ e -1/2 3/2 

Y B -l(oo) a 1 z 



Zr exp 
9* 



-Zf — / dzP(z) 

J Zf 



(26) 



where a = —dP(z)/dz\ z=Zf (> 0) with Zf defined as P(zf) = 1. In the first case, -Pd(I) 3> 1 
but Pjv(1) "C 1, it is difficult to obtain a simple formula because of the Bessel function 
in Pd(z). Although we can find a rough formula Yb-l cc 1/k, more reasonable fitting 
formula is known such as jn| 

Y °-^ ^ i x Mi^r (27) 

with k ~ -Pd(I)- In the other case, -Pd(I) 1 but Pjv(1) 3> 1, analytic calculation is 
possible by using Eq. (pl~D, and the approximation formula is found to be 

Yb-l(oo) ~ — x Aexp (-2v / A) (28) 

with 



47r 3 F(M 1 )^ 



(29) 



We can check that the both formulas give good approximations compared to the numerical 
results. 
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We have learned that there is the parameter region where the washing-out processes 
become important. Here, remember that e is the function of the Yukawa couplings and 
becomes large according to the power law, e oc h 2 . In the case -Pd(I) 1 but .P/v(l) <C 1, 
from Eq. (|27|), we can expect that the resultant baryon number is saturated to a con- 
stant, as Yukawa coupling becomes large. On the other hand, in the case -P/v(l) ^> 1 but 
-Pd(I) ^ 1) Eq. ( p8|) leads to the conclusion that the resultant baryon number becomes 
an exponentially decreasing function of the (large) Yukawa couplings. This implies that, 
even if the parameter e/g* itself is too large to be consistent with the observation, new pa- 
rameter region of large Yukawa couplings comes out, so that the observed baryon number 
can be reproduced. This is a new interesting scenario. 

In summary, we have studied the leptogenesis scenario in models with multi-Higgs 
doublets. In our scenario, the direct correspondence between the Yukawa coupling and 
the mass matrix is lost. As a result, new parameter region can be consistent with the 
neutrino oscillation data. We have pointed out the importance of the washing-out process 
which has not been seriously taken into account in the conventional scenario. We found 
a new parameter region of the large Yukawa couplings which can reproduce the observed 
baryon number through the exponential dumping of too large CP-violating parameter e 
due to the washing-out processes through the dimension five interactions . 

Finally, we give a comment on the application of our scenario to some neutrino mass 
matrix models. Many models of the Majorana neutrinos with the see-saw mechanism 
have been proposed, which predict and/or reproduce the observed neutrino properties. 
When the leptogenesis scenario is taken into account in such models, some of them may 
predict too large baryon asymmetry to be consistent with the observations. This is very 
problematic in cosmology.^ However, if such a model is the multi-Higgs doublet model, 
the scenario we have discussed in this letter is applicable and can make the model viable. 
Further discussion needs a concrete model, and such a model is worth investigating. 

The authors would like to thank Rabindra Mohapatra and Markus Luty for useful 
discussions. 



4 If we assume that the re-heating temperature of the inflationary universe is lower than the right- 
handed neutrino masses, this problem can be easily avoided. However, the leptogenesis no longer works. 
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Figure 1: The solution of the Boltzmann equation (Log 10 (Y B _ L x 10 10 )) with e/g* = 
10~ 5 (the upper horizontal line), M\ = 10 10 GeV and rn u3 = 0.2 eV. The solutions for 
m ul = 4 x 10~ 3 , 0.04, 0.4, 4 and 40 eV are plotted from above at z = 10 2 . 
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Figure 2: The solution of the Boltzmann equation (Log 10 (Ys_L x 10 10 )) with e/ g* = 10~ 5 , 
Mi = 10 10 GeV and m ul = 4 x 10~ 4 eV. The solutions for the input values fn u3 = 7, 12, 
17, 22 and 27 eV are plotted from above. 
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